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Scalar and spin-dependent relativistic effects on magnetic properties calculated with 
four-component methods: the nuclear magnetic resonance parameters of the lead 

halides 
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The results of calculations of nuclear magnetic resonance (NMR) parameters for the lead halides 
is reported in this paper. The results are obtained by using four-component methods. The use of 
the nonrelativistic Levy-Leblond Hamiltonian along with the relativistic Dirac-Coulomb and spin- 
free ones allows us to discriminate scalar and spin-dependent effects on the parameters. It is found 
that the wide range of the lead NMR spectra and their large anisotropies are, mainly, due to 
spin-dependent effects on the paramagnetic term. Among the relativistic scalar corrections, the so- 
called spin-Zeeman kinetic-energy term turns out to be dominant. The reduced spin-spin coupling 
constants become proportional to the product of the atomic numbers of the coupled nuclei. 
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I. INTRODUCTION 



It is nowadays well established that the calculation of 
properties of molecules containing heavy atoms requires 
the inclusion of the kinematic effects of relativity in or- 
der to correctly account for their magnitude and, in some 
cases, even their experimental trends [ll 3- A number 
of methods and formalisms have been devised for such a 
purpose 1^, [3!, \4, l5|, IS!, |7|, |8| ; among them, a6 initio fully 
relativistic four-component methods [1, 0, [a] occupy a 
privileged position as a reference for benchmark calcu- 
lations against which more approximated methods are 
to be compared. In particular, the calculation of the pa- 
rameters of the nuclear magnetic resonance (NMR) spec- 
troscopy are known to be highly sensitive to the relativis- 
tic effects 9, 10]. 

From an experimental point of view, NMR spec- 
troscopy is a powerful tool for the determination of 
molecular structures, thus enabling its use for a wide 
range of applications. Among the magnetic nuclei, ^°^Pb, 
the heaviest one having spin / — 1/2, is particularly suit- 
able for discriminating among different structures and ox- 
idation states in lead-containing compounds because of 
its natural abundance (22.6%) and excellent receptivity 
[9|| . The high sensitivity of ^''^Pb to its coordination envi- 
ronment gives rise to large changes in its chemical shifts 
covering a range of about 16000 ppm for the known com- 
pounds. Lead can form both inorganic and organolead 
compounds. Inorganic lead chemistry is dominated by 
divalent Pb(II) while organolead chemistry is based on 
tetravalent Pb(IV) [ll|. Lead dihalides are among the 
simpler inorganic salts that divalent lead can form. PbF2 
is probably the most studied lead halide. Solid PbF2 
is a ionic solid presenting a strong superionic conduc- 
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tivity. It exhibes two phases: an orthorombic a-PbF2 
phase stable at high pressure and low temperatures, and 
a cubic /9-PbF2 phase more stable in ambient conditions. 
[l3 |. The absoption spectra of gaseous PbF2 was first 
reported in Ref. [13|. It was shown that molten PbF2 
vaporizes partly as molecular PbF2 with a considerable 
desproportionation jlJI. Measurements of ^°^Pb chem- 
ical shifts in lead halide powders have been performed 
using the magic-angle spinning (MAS) technique Il4|. 
Calculations of ^^F NMR shielding constants has been 
reported for a number of metal fluorides, including PbF2 
[la | . Calculations of molecular structures and vibrational 
spectroscopic properties for the lead halides, among other 
group 14 compounds, have been performed by using 
a correlated complete active space self-consistent field 
(CASSCF) and multireference singles-doubles configura- 
tion interaction (MRSDCI) wave function, and relativis- 
tic effective core potentials (ECP) [la, U^, the coupled 
cluster method with singles and doubles excitations and 
perturbative inclusion of triple excitations, CCSD(T), su- 
plemented with ECPs [l8| , and density functional theory 
(DFT) and ECPs ^. To the author's knowledge, the 
only theoretical study of the NMR magnetic shieldings 
of isolated molecules of lead halides is the one presented 
in Ref. [2^ . In that paper it is also reported calculations 
in clusters that preserves the orthorombic symmetry of 
the solid and experimental measurements in solid sam- 
ples [20|. However, although an agreement is found with 
experimental trends, the basis sets used are small and the 
method of calculation involves too many approximations 
so that no precise quantitative conclusions can be drawn 
from the results. 

In this paper we shall consider the fully relativistic 
calculation of the NMR parameters of the lead dihalides 
with a twofold purpose. Firstly, it is aimed to explain 
the wide-range variations of the NMR parameters, their 
large anisotropies, and their experimentally observed reg- 
ularities as the different halogens are substituted, as well 
as to identify the underlying mechanisms beneath such a 



behavior. Secondly, the use of all-electron fully relativis- 
tic methods along with large-size basis sets, allows us to 
test the performance of these reference methods in sys- 
tems of analogous molecules ranging from PbF2 , in which 
the central heavy atom is dominant, to Pbl2, constituted 
by three heavy atoms. 



II. CALCULATION APPROACH 
A. Computational details 

Fully relativistic, spin- free and nonrelativistic calcu- 
lations of the electronic structure of the lead dihalides 
PbX2 {X = F, CI, Br and I) were carried out self- 
consistently with the Dirac-Coulomb, spin-free [2l| and 
Levy-Leblond Hamiltonians |22l |. respectively. The nu- 
clear magnetic shielding tensor, a, and the reduced 
spin-spin coupling tensor K were calculated as four- 
component linear response functions using the polar- 
ization propagator approach within the Random Phase 
Approximation (RPA) Q. All calculations were per- 
formed with the DIRAC code [23 . It should be noted that 
even the nonrelativistic calculations were done with four- 
component methods by using the nonrelativistic Levy- 
Leblond Hamiltonian. Usually, fully relativistic and non- 
relativistic calculations require different basis sets be- 
cause of the different nature of the relativistic and non- 
relativistic operators. All along the calculations reported 
here, the relativistic bases were used and have shown 
to provide convergence both in the spin-free and Levy- 
Leblond calculations. The use of a large value of the 
speed of light in fully relativistic calculations has also 
been advocated as a procedure to reach the nonrelativis- 
tic limit. 

The molecular structures were taken from experimen- 
tal geometries [HI,!!!: R{X-Vh) = 2.033, 2.46, 2.6 and 
2.79 A, and e{X-Vh-X) = 97.8, 96.0, 98.8 and 99.7° for 
X = F, CI, Br and I, respectively. These bond lengths 
are close to values calculated with the MRSDCI method 
[la |. but the calculated bond angles are about four de- 
grees larger than the ones used here. So, there is some 
uncertainty inherent to our calculations because of their 
dependence on the molecular geometry. In this work, 
large component uncontracted Cartesian Gaussian basis 
set were used in all atoms. The large component basis 
sets used in the calculations were obtained from modi- 
fying the Faegri basis sets (for Pb, Br and I) [2^ and 
the aug-ccpVTZ basis set (for F and CI) by suplement- 
ing them with large- (tight functions) and small-exponent 
(diffuse functions) Gaussian basis functions. These tight 
(diffuse) functions were obtained from multiplying (di- 
viding) by a factor of 3 the largest (smallest) exponent 
present in each block of angular momentum. The small 
component bases sets are generated from the large com- 
ponent bases by invoking the kinetic balance condition. 
The final ranges of the exponents in the expanded bases 
are similar to those used previously in accurate nearly 



basis set-limit calculations [26 



B. Basis set optimization in PbBr2. 

Given that the basis set completeness is a critical is- 
sue for reliable results, specially in the calculations of the 
two heaviest halides, the PbBr2 molecule was taken as a 
model to test the basis set quality as a function of the 
functions added to the original basis sets for the large 
component. Table [J shows the calculated values for the 
paramagnetic, diamagnetic and total nuclear magnetic 
shielding in both the lead and bromine nuclei and the 
one- and two-bonds indirect spin-spin coupling constants 
V(20^Pb-™Br) and ^(^^Br-^Br). As one goes from ba- 
sis A (Faegri basis set) to the large-size basis G, all the 
quantities converges steadily. In particular, the shielding 
diamagnetic terms a'^ (Pb) and a'^ (Br) and the two-bonds 
coupling ^J(™Br-^^Br) are reasonably well reproduced 
with the basis A already. In contrast, the shielding para- 
magnetic terms and the one-bond coupling ^J(^"^Pb- 
^^Br) require of the large bases F or G to achieve conver- 
gence within about 10 ppm for aP(Pb), less than 1 ppm 
for crP(Br) and ca. 3 Hz for ^ J(Pb-Br). Since no attempt 
was done to include electronic correlation effects, the ba- 
sis F was considered converged enough and it was used 
for representing Pb and Br in all the calculations below. 
A similar basis scheme was adopted for the other halo- 
gens, including functions with exponents of the same or- 
der of magnitude as those used in similar atomic calcula- 
tions [2^. The effect of the tight, diffuse and polarization 
functions added was studied, although not shown here, 
to ensure that convergence was reached. With such cri- 
teria, the PbF2, PbCl2, PbBr2 and Pbl2 molecules were 
described by 2168, 2288, 2332 and 2343 Gaussian basis 
functions, respectively. The large and small component 
basis sets finally used are listed in Table [TTl 



C. Neglect of (S'5'|5'S)-integrals 

Following an approximation proposed by Visscher [27| | , 
the contribution from (S'5'|S'S')-type two-electron inte- 
grals, i.e. those involving small components functions 
only, were neglected both in the wave function and the 
response calculations along this paper. This approx- 
imation allows a considerable saving in computational 
demands, typically decreasing the timings in one order 
of magnitud. However, an estimation of their impor- 
tance was obtained by calculating their contributions 
in the PbBr2 molecule with the smallest basis set A. 
The use of this basis set instead of the larger ones is 
justified because the small components of the molecu- 
lar orbitals are densely centered in the atoms and only 
affect the electronic density in the very proximity of 
the nuclei, in contrast to the diffuse functions added 
to basis set A in the bases from B to G. Because of 
that, quite different combination of basis sets [23, l2g| 



TABLE I: Calculated values of the NMR parameters for PbBr2 as a function of the basis set. 



Large component criso(Pb) criso(Br) 

Basis functions No. param. diam. total param. diam. total ^J(Pb-Br) ^J(Br-Br) 



Pb: 23s 21p Ud 9/ 
A Br: 19s 16p 9d 

Pb: 27s 22p 17d 12/ 2g Ih 
B Br: 19s 16p 9d 1/ 1855 -8119.8 8025.1 



1511 -9466.4 7974.2 -1492.2 -216.2 2606.3 2390.1 

-94.7 -498.5 2624.9 2126.4 



Pb: 27s 22p 17d 12/ 2g Ih 
C Br: 19s 16p 9d 2/ 1905 -8190.0 8026.1 -163.1 -529.8 2628.6 2098.8 

Pb: 27s 22p 17d 12/ 2g Ih 
D Br: 19s 16p 9d 3/ 1955 -8196.5 8027.4 -169.9 -530.3 2630.7 2100.4 

Pb: 27s 22p 17d 12/ 2g Ih 
E Br: 20s 16p 13d 2/ Ig 2125 -8223.4 8028.0 -195.3 -545.3 2668.5 2123.2 



Pb: 27s 22p 17d 12/ 3^ 2h 

Br: 20s 16p 13d 3/ 2g 2332 -8122.9 8030.9 

Pb: 27s 22p 17d 12/ 4g 3h 

Br: 20s 16p 13d 3/ 2g 2417 -8113.9 8031.1 



-92.0 -551.1 2670.3 2119.2 
-82.7 -551.8 2670.6 2118.8 



-3446.63 


43.26 


-3847.67 


46.77 


-3866.52 


46.18 


-3868.19 


46.13 


-3891.24 


46.16 


-3950.06 


44.99 


-3952.79 


44.94 



like Faegri(Pb)/Faegri(Br), Faegri(Pb)/Sadlej(Br) and 
Hirao(Pb)/Sadlej(Br) yield quite different i7iso(Pb) and 
criso(Br) far away from the converged values but, never- 
theless, all of them give almost the same contributions 
from the inclusion of (S'S'|55)-type two-electron inte- 
grals, namely, ~-|-50 ppm to CTiso(Pb) and ca. +1 ppm 
to criso(Br). Hence, their contributions to (Tiso(F) and 
criso(Cl) are negligible, while they contribute a few ppm 
to criso(I)- All these contributions are below 1% of the 
total values, thus justifying the exclusion of these in- 
tegrals. Another approximation consisting in including 
the (iL|S'S')-type integrals in the wave function calcula- 



J 



tion although neglecting them in the response calculation 
of the properties was also tried with similar results. It 
turns out that they typically contribute some ~ +40 ppm 
to criso(Pb), somewhat less (~ -1-10 ppm) for criso(I) and 
their effect is negligible for the lighter halogens. This 
approximation yields an extra saving in computing time 
and it could be recommendable in order to speed up com- 
putations without too much lost of accuracy. However, 
for consistency, all the numbers reported in this paper 
were obtained by fully including {LL\SS) two-electron 
integrals in every step of the calculation. 
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TABLE IL Large and small component basis sets used in the 
calculations. 



Atom 


Large component 


Small component 


F 


13s lip 9d 6/ 


lis 22p 17d 9/ 63 


CI 


20s 13p 8d 7/ 


13s 28p 20d 8/ 7g 


Br 


20s 16p 13d 3/ 2g 


16s 22p 16d 13/ ig 2h 


I 


23s 20p 16d 2/ Ig 


20s 28p 20d 16/ 2g Ih 


Pb 


27s 22p 17d 12/ 3g 2h 


22s 29p 25d 17/ 12^ 3/i 2i 



III. RESULTS AND DISCUSSION 
A. Absolute nuclear shieldings 

The results of nuclear magnetic shielding tensors for 
lead and halogen nuclei in the series of molecules PbF2, 
PbCl2, PbBr2 and Pbl2 calculated with relativistic, spin- 
free and nonrelativistic Hamiltonians are shown in Table 
IIHI The tensor cr{Fh) is diagonal, with components fTii, 



a±_ and Uz along the direction parallel to the molecular 
symmetry axis, perpendicular to it but in the plane of 
the molecule, and perpendicular to the molecular plane, 
respectively. In Table IIHI it is also shown the individual 
components of the paramagnetic and diamagnetic tensors 
crP and <t'^, their corresponding isotropic values tTiso = 
(l/3)Tr(<T), and the span fi == ai — 0-3 defined as the 
diference between the largest and the smallest eigenvalues 
of the tensor (cti > 172 > cr^) for Pb and the halogen 
atoms. 

Table IIVI provides further insight on the results shown 
in Table IIHI by decomposing the various magnitudes in 
scalar and spin-dependent contributions given by 

(1) 

_spin — free /ir\\ 



calar 



SO 



.^spin— free .^NR 



^Rel 



such that a = a 



^NR 



^scalar 



a^'~' . A similar decom- 



position for the spin-dependent terms has been recently 
applied to study of the static polarizability of group- 13 
atoms |29|]. 



Pb 



A clear trend to deshielding {a < 0) of the lead nu- 
cleus in PbX2 is observed from the relativistic (Rel), the 
spin-free (SF) and the nonrelativistic (NR) calculations 
as one goes from X =F to CI to Br to I. However, this 
trend is much more pronounced for tTp'^'(Pb) than for 
tTj^Q (Pb) or (7?^^(Pb). The analysis of aiso in terms of 
its paramagnetic and diamagnetic contributions shows 
that this is a consequence of the dramatic increase of the 
(negative) paramagnetic term from ^ -f 200 (in PbF2) to 
^ —14000 ppm (in Pbl2), that cannot be compensated 
by the slight growth in af^^ from 7896 to 8110 ppm. Fur- 
thermore, an inspection of the tensor components of the 
relativistic ctP and cr'^ shows that, to a large extent, cr"^ 
is almost isotropic, reflecting its deep core origin, while 
the different components of the tensor <tP behave dis- 
tinctly in response to the kinematic effects of relativity, 
CT^. being more sensitive than ctP, and this, at its turn, 



more sensitive than a?^, yielding a^ < a^ < a 



II 



all 



cases, both relativistically and nonrelativistically. It has 
been noted that the measured ^°^Pb isotropic shifts in the 
lead dihalides fit to a linear correlation with the inverse 
of the ionization potential fl4|; such a correlation has 
been taken as an indication that the paramagnetic con- 
tribution to the chemical shift is dominant. Reciprocally, 
given that in our calculations the paramagnetic contribu- 
tions arc the most affected by relativistic effects, and tak- 
ing the HOMO-LUMO gap as a measure of the ionization 
potential, Figure [T] shows that there is a linear correlation 
between each tensor component of the relativistically cal- 
culated crP(Pb) with the inverse of enoMO-LUMO- The 
different slopes of the least-squares fitted straight lines 



can be taken as proportional to the corresponding ma- 
trix elements of the perturbations giving a?^, a^ and a^. 



2. halogen atom shieldings 



As shown in Table|lVl both a^^^^ {X) and a] 



p/scalar 



{X) 



becomes more positive as one changes the sustituent from 
F to I. The relative importance of their magnitude, how- 
ever, changes from I^Pf °(X)| > W\^^''''''\X)\ for X =F 
and CI, to the opposite for Br and I. There is no obvi- 
ous interpretation for such a trend. In previous works, 
the study of the correlation of the various contributions 
to the nuclear charge Z has been a useful guide for the 
identification of underlying mechanisms. The best least- 
squares fits found for the paramagnetic shielding con- 
stants at the halogen nuclei is given by: 



C: 



p/ scalar 
p/SO 



{X) 
{X) 



3.19 X IQ-'^Z^"^^ - 


-27.8 


(3) 


1.12x10-5^4-29- 


-41.3 


(4) 



Two facts deserve some attention: firstly, the different 



Z-dependence of those contributions, i.e., CT: 



p/scala 



Z3 



while a] 



p/SO 



Z''; and secondly, these correlations of the 



form {AZ^ + B) fits better to the calculated values than 
the simpler ones {AZ") which one would expect from a 
relativistic correction, in order to fulfill the requierement 
of vanishing when Z goes to zero. No further insight 
can be provided within the present framework. A con- 
sideration of a perturbative treatment could be helpful 
in analyzing these trends more deeply. 



r 



'Pb shifts 



The standard reference for ^^^Pb chemical shifts mea- 
surements is tetramethyllead (TML), Pb(CH3)4, i.e., 
the chemical shifts are defined as S{'^°''Ph)= cr(TML) — 
cr(PbAr2). Nevertheless, we shall estimate here the lead 
chemical shifts by refering our compounds to plumbane 
(PbH4), i.e., with respect to a reference in which every 
methyl group was replaced by a hydrogen atom; such a 
change should not introduce severe changes while keep- 
ing the computational costs withing reasonably bounds. 
The calculation of plumbane NMR parameters is not only 
useful for obtaining the measurable shifts but also be- 
cause of its tetrahedral symmetry T^, what provides an 
interesting example for the analysis of the spin-dependent 
and scalar contributions to the nuclear magnetic shield- 



ing tensor of lead. Because of the Td symmetry, the 
tensor cr (Pb) as well as its paramagnetic and diamag- 
netic contributions become isotropic. The results of 
Dirac-Coulomb, spin-free and nonrelativistic calculations 
are shown in Table |Vl It can be seen that the spin- 
dependent contributions to both the paramagnetic and 
diamagnetic term represent a small fraction of the to- 
tal value, while almost the total relativistic corrections 
are due to scalar effects. Furthermore, the magnitude of 
those effects (crfi'"''''"' c^ 8000 ppm, aH^"''^''' ~ -2300 
ppm) are quite close to, e. g., those in PbF2, thus show- 
ing that they arise from the inner most electrons in the 
core. The isotropic lead chemical shifts and its param- 
agnetic and diamagnetic contributions are also shown in 
Tabic |Vl The decomposition of the relativistic correc- 
tions to the shifts {5^°^ - 5^^) in terms of scalar (J'^'^aiar) 



TABLE III: Isotropic values and tensor components of the nuclear magnetic shielding tensor of the lead and halogens nuclei X 
{X = F, CI, Br, I) calculated with relativistic (Rel.), spin-free (SF) and nonrelativistic Levy-Leblond (NR) Hamiltonians. All 
values are given in ppm. 



Nucleus 


Quantity 


Rel. 


PbF2 
SF 


NR 




PbCl2 






PbBr2 






Pbl2 






Rel. 


SF 


NR 


Rel. 


SF 


NR 


Rel. 


SF 


NR 


^"''Pb 


CTiso 


8078.5 


13226.5 


7465.4 


3152.0 


11317.4 


6101.3 


-92.0 


10511.7 


5616.1 


-5970.4 


9273.7 


4763.8 




<o 


182.2 


5371.6 


-2676.3 


-4790.8 


3416.0 


-4086.7 


-8122.9 


2519.7 


-4695.0 


-14080.7 


1202.8 


-5651.4 




-L 


7896.4 


7854.9 


10141.7 


7942.9 


7901.4 


10188.1 


8030.9 


7992.0 


10311.1 


8110.2 


8070.9 


10415.2 




o"! 


12723.9 


14233.1 


8561.4 


9120.4 


12463.7 


7463.3 


6239.6 


11521.9 


6954.7 


653.1 


10105.2 


6125.8 




CTz 


5974.1 


12721.5 


7083.4 


2427.1 


11420.7 


6216.2 


559.2 


10982.3 


6078.7 


-1944.0 


10640.8 


6041.3 




0-x 


5537.5 


12725.0 


6751.4 


-2091.4 


10067.8 


4624.5 


-6974.9 


9030.8 


3815.0 


-16620.4 


7075.2 


2123.8 




n 


7186.4 


1508.1 


1810.0 


11211.8 


2395.9 


2838.8 


13214.5 


2491.1 


3139.7 


17273.5 


3030.0 


4002.0 




''\\ 


4840.5 


6391.3 


-1567.0 


1199.5 


4584.5 


-2702.4 


-1762.2 


3558.8 


-3322.3 


-7421.3 


2070.1 


-4245.6 




ol 


-1963.7 


4825.0 


-3100.0 


-5578.1 


3456.8 


-4034.7 


-7676.8 


2884.6 


-4355.4 


-10196.7 


2426.2 


-4546.7 




< 


-2330.3 


4898.6 


-3361.8 


-9994.0 


2206.7 


-5523.1 


-14929.7 


1115.7 


-6407.3 


-24624.0 


-887.8 


-8161.8 




<rt 


7937.9 


7896.5 


10183.4 


8005.2 


7963.9 


10250.9 


8136.0 


8097.6 


10434.1 


8252.7 


8214.6 


10588.3 






7883.4 


7841.8 


10128.4 


7920.9 


7879.2 


10165.7 


8001.8 


7963.1 


10277.0 


8074.4 


8035.1 


10371.4 




7867.8 


7826.3 


10113.2 


7902.6 


7861.1 


10147.6 


7954.8 


7915.2 


10222.3 


8003.6 


7963.0 


10285.7 


X 


O"iso 


220.3 


263.1 


289.0 


641.4 


678.6 


750.5 


2119.2 


2113.7 


2260.1 


4453.7 


4210.4 


3980.6 




-L 


-261.7 


-219.0 


-197.9 


-448.2 


-410.8 


-415.2 


-551.1 


-556.2 


-904.3 


91.8 


-142.4 


-1580.4 




-L 


482.0 


482.1 


486.9 


1089.6 


1089.4 


1165.7 


2670.3 


2669.8 


3164.4 


4361.9 


4352.8 


5561.0 




o\ 


458.9 


462.0 


463.7 


994.1 


993.2 


1023.6 


2812.4 


2810.1 


2831.6 


5397.1 


5357.5 


4812.6 




02 = 0-^ 


219.3 


268.0 


299.5 


766.4 


794.9 


874.7 


2494.7 


2452.0 


2613.7 


5355.9 


4991.6 


4750.4 




0-3 


-17.2 


59.4 


103.9 


163.6 


246.6 


353.3 


1050.5 


1078.9 


1335.1 


2608.1 


2282.1 


2378.7 




n 


476.1 


402.6 


359.8 


830.5 


746.6 


670.3 


1761.9 


1731.2 


1496.5 


2789.0 


3075.4 


2433.9 




'y\ 


-50.6 


-48.5 


-53.0 


-118.6 


-120.5 


-168.8 


116.3 


115.6 


-364.9 


1008.7 


980.6 


-784.7 




Tl 


-258.8 


-210.3 


-183.0 


-327.2 


-298.7 


-294.6 


-187.7 


-230.1 


-565.3 


974.2 


618.6 


-834.1 




<^l 


-475.5 


-398.2 


-357.6 


-899.0 


-815.0 


-782.3 


-1581.9 


-1554.0 


-1782.8 


-1707.4 


-2026.4 


-3122.4 




at 


509.6 


509.7 


515.5 


1112.7 


1112.7 


1189.7 


2696.3 


2696.1 


3194.6 


4390.7 


4382.1 


5597.2 




at 


478.1 


478.2 


482.5 


1093.6 


1093.6 


1169.2 


2682.4 


2682.1 


3178.9 


4381.6 


4373.0 


5584.5 




ai 


458.3 


458.5 


462.7 


1062.6 


1062.6 


1138.2 


2632.2 


2631.3 


3119.8 


4313.2 


4303.2 


5501.2 



and spin-dependent term J^^ shows that almost the en- 
tire correction is due to the spin-orbit contribution to the 
paramagnetic term iJp/^O in contrast to (5'^/^'-' which be- 
comes negligible because 6'^ comes mainly from contribu- 
tions from the spherically symmetric inner most atomic- 
hke s orbitals. The magnitude of the contributions be- 
comes (5^^/^° < ^d/scalar ^ Jp/scalar ^ ^p/SO^ ^^^^ ^p/SO 

accounting for more than 90% of Siso even for Pbl2. 



C. Reduced spin-spin couplings 

The inclusion of relativistic effects in the calculation of 
indirect spin-spin coupling constants of lead-containing 
compounds has been considered through a number of 
methods, e.g., semiempirically [30, l3l|j DFT-ZORA ap- 
proximation |32l . |33| and fully relativistic polarization 
propagator approach [SJI, although, to the author's 
knowledge, there are no calculations for the couplings 
in the lead dihalides. In Table FvH the results of nonrela- 
tivistic (NR) and relativistic (Rel) calculations are shown 
along with the scalar and spin-dependent (SO) contribu- 
tions. As with the a and S calculations, K^'^^^'^'^ and K^'^ 



are obtained as a difference with respect to the result of 
a spin-free calculation. It can be seen that both isotropic 
and anisotropic one-bond couplings becomes more and 
more negative as one goes from the fuoride to the iodide. 
Both the scalar and spin-dependent contributions give 
rise to such a trend, with the scalar mechanisms being 
dominant. 

On the contrary, two-bonds couplings show a decreas- 
ing behavior in nonrelativistic calculations, as opposite 
to the increasing trend of the relativistic ones, as one 
goes from F to CI to Br to I. This is a consequence of, 
mainly, the SO and -to less extent- scalar mechanisms, 
both of which give positive contributions to ^i^iso and 
A^K. Nevertheless, it is worthy to note that, for the 
heaviest halide, Pbl2, A^ K^'^'^^'^" (l—V) becomes greater 
than A^ K^'-' (I—T) . There have been attempts in the lit- 
erature to correlate the reduced spin-spin coupling con- 
stants with the square of the atomic number of a central 
heavy atom bonded to light atoms. On the other hand, 
Bryce et al. found that Kiso and AK have a good linear 
correlation with the product of atomic numbers of the 
coupled nuclei in interhalogen diatomics [32l | . In the Fig- 
ure[2]we plot relativistically calculated ^i^iso and A^K as 



TABLE IV: Isotropic values and tensor components of the spin-orbit and scalar contributions to the nuclear magnetic shielding 
tensor of the lead and halogens nuclei. All values are given in ppm. 



Nucleus Quantity 


PbF 


'2 


PbCl 


2 


PbB: 


r2 


Pbl2 






SO 


Scalar 


SO 


Scalar 


SO 


Scalar 


SO 


Scalar 


""'Pb aiso 


-5148.0 


5761.1 


-8165.4 


5216.1 


-10603.7 


4895.6 


-15244.1 


4509.9 


<o 


-5189.4 


8047.9 


-8206.8 


7502.7 


-10642.6 


7214.7 


-15283.5 


6854.2 




41.5 


-2286.8 


41.5 


-2286.7 


38.9 


-2319.1 


39.3 


-2344.3 


'^ll 


-1509.2 


5671.7 


-3343.3 


5000.4 


-5282.3 


4567.2 


-9452.1 


3979.4 


(Tz 


-6747.4 


5638.1 


-8993.6 


5204.5 


-10423.1 


4903.6 


-12584.8 


4599.5 




-7187.5 


5973.6 


-12159.2 


5443.3 


-16005.7 


5215.8 


-23695.6 


4951.4 


"l 


-1550.8 


7958.3 


-3385.0 


7286.9 


-5321.0 


6881.1 


-9491.4 


6315.7 


<T? 


-6788.7 


7925.0 


-9034.9 


7491.5 


-10561.4 


7240.0 


-12622.9 


6972.9 


< 


-7228.9 


8260.4 


-12200.7 


7729.8 


-16045.4 


7523.0 


-23736.3 


7274.0 


<rt 


41.4 


-2286.9 


41.3 


-2287.0 


38.4 


-2336.5 


38.1 


-2373.7 


1 
(7± 


41.6 


-2286.6 


41.7 


-2286.5 


38.7 


-2313.9 


39.3 


-2336.3 


41.5 


-2286.9 


41.6 


-2286.6 


39.6 


-2307.1 


40.6 


-2322.7 


X (Tiso 


-42.8 


-25.9 


-37.2 


-71.9 


5.5 


-146.4 


243.3 


229.8 


<so 


-42.8 


-21.1 


-37.4 


4.4 


5.0 


348.1 


234.2 


1438.0 




-0.1 


-4.8 


0.2 


-76.3 


0.5 


-494.6 


9.1 


-1208.2 


cri 


-3.1 


-1.7 


0.9 


-30.4 


2.3 


-21.5 


39.6 


544.9 


02 = 0-2 


-48.7 


-31.5 


-28.5 


-79.8 


42.7 


-161.7 


364.3 


241.2 


0-3 


-76.6 


-44.5 


-83.0 


-106.7 


-28.4 


-256.2 


326.0 


-96.6 


n 


73.5 


42.8 


83.9 


76.3 


30.7 


234.7 


-286.4 


641.5 


'^f 


-2.1 


4.5 


1.9 


48.3 


0.7 


462.5 


28.1 


1765.3 


<rl 


-48.5 


-27.3 


-28.5 


-4.1 


42.4 


335.2 


355.6 


1452.7 


'^l 


-77.3 


-40.6 


-84.0 


-32.7 


-27.9 


228.8 


319.0 


1096.0 


<rt 


-0.1 


-5.8 


0.0 


-77.0 


0.2 


-498.5 


8.6 


-1215.1 


T^ 


-0.1 


-4.3 


0.0 


-75.6 


0.3 


-496.8 


8.6 


-1211.5 


4 


-0.2 


-4.2 


0.0 


-75.6 


0.9 


-488.5 


10.0 


-1198.0 



a function of the halogen atomic number Zx (left panel) 
and ^i^iso and A^K as a function of Z\ (right panel 
panel). The dashed line represents the least-squares fit 
of the data. As it can be seen, the correlation is quite rea- 



J 



sonable, especially for the isotropic values and somewhat 
less for the anisotropics. These results provide support 
to the correlation proposed in Ref. [32| . generalizing it 
to include two-bonds couplings. 



r 



IV. CONCLUSIONS 

In summary, in this paper we have reported results of 
calculations of the NMR parameters for the lead halides. 
The calculations made use of four-component relativistic 
Dirac-Coulomb and nonrelativistic Levy-Leblond Hamil- 
tonians. In addition, the use of the spin-free hamiltonian 
allows to separate the scalar and spin-dependent rela- 
tivistic contributions. 

The results presented in this paper can be rationalized 
as follows: 

(i) the main relativistic correction to the nuclear mag- 
netic constants in lead is the spin-dependent contribution 
(jP/SO j-pj^-j which can be roughly associated with the spin- 
orbit interaction, strongly affected by the valence struc- 
ture. Substitution of increasingly heavy halogens changes 



dramatically such a structure from PbF2 dominated by 
the lead nucleus to Pbl2, a three heavy atoms molecule. 

(ii) o-P/'^'=^i^'(Pb) ranges from 8047 ppm for PbF2 to 
6854 ppm for Pbl2, what is compatible with a dominant 
^sz-KE |;(;pj-Q estimated in 7414 ppm by using an ana- 
lytical model for the shielding in atoms [35|- However, 
the corresponding estimations for F (9 ppm) and CI (64 
ppm) are not good because of the large anisotropics of 
their o-p/scaiar tensors; the model estimates 570 and 1995 
ppm for Br and I, in moderate agreement with the val- 
ues of Table|lV](348 and 1438 ppm). The reason for that 
is that the spherical symmetry assumed in the model is 
more approximately realized in the heaviest halogen nu- 
clei. 

(in) The (j'^h'^^^^'^ term was studied in Ref. [3^| , where 
it has been shown that, within a perturbative frame- 
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FIG. 1: Isotropic and tensor components of the paramagnetic contribution to the relativistic Pb absolute nuclear magnetic 
shielding in (from left to right) PbF2, PbCl2, PbBr2 and Pbl2, as a function of the inverse of the HOMO-LUMO gap. 



work, it arises from two linear response contributions 
(Darwin and mass-velocity) plus three expectation value 
terms, with all of them following a Z^ power-law scal- 
ing with the atomic number; e.g.^ the perturbative LR- 
ESC method gives a'^/^''''^^'{Rn) = -3446.4 ppm ^, 
whose extrapolation to Pb would give cr'^/'''^^''''^(Pb) = 
-3446.4 X (82/86)3 = -2987.5 ppm, in line with the 
value of about -2300 ppm obtained for all halides. Per- 
forming the same extrapolation from the LR-ESC value 
for the nearest rare gas atom (taken from Ref. ^0|) to 
each halogen, one has -3.3, -23.8, -215.1 and -767.1 ppm, 
to be compared to -4.8, -76.3, -494.6 and -1208.2 ppm 
of the present work (Table ITV)) . Although the trend is 
correct, a lack of quantitative agreement still remains. 

(iv) The cr'^/^*-' term can be disregarded in all cases, 
what can be attributed to the negligible effect of the 
spin-orbit interaction upon the inner electronic shells, 
which are responsible for the diamagnetic shielding. 

The large spin-dependent effects on the paramagnetic 
shielding are also responsible for the large observed vari- 
ations of the lead resonance over a wide range, as well 
as its large anisotropics, which originates in the different 
responses of the molecule to the spin-orbit interactions, 
although every component become proportional to the 



inverse of the HOMO-LUMO gap. The same trends 
are found in the ^"^Pb chemical shifts in the various 
halides because of the isotropic nature of the relativistic 
corrections, having a scalar origin, in the reference 
compound (plumbane). 

Both the scalar and the spin-dependent mechanisms 
are important for the isotropic reduced coupling con- 



stants ^Kiso and ^K\ 
while pifgjf'-l < , .-,,, 



.,o, although |ii^£Jf'-'-| > |iif?0|, 
^KpJ^\. Finally, the relativistic one- 



and two-bonds reduced coupling constants becomes well 
correlated to the product of the nuclear charges of the 
coupled nuclei, Zp\^Zx and Zj^, respectively. 
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TABLE V: Isotropic nuclear magnetic shielding constant of 
lead nucleus in plumbane calculated with relativistic (Rel.), 
spin-free (SF) and nonrelativistic Levy-Leblond (NR) Hamil- 
tonians, and spin-orbit (SO) and scalar contributions to them, 
and isotropic chemical shifts of the Pb nuclei in the lead 
halides with respect to plumbane, taken as reference com- 
pound. All values are given in ppm. 



Molecule 




Rel. 


SF 


NR 


SO 


scalar 


PbH4 


a 


12776.7 


12674.2 


6869.5 


102.5 


5804.7 




o-P 


4941.3 


4880.4 


-3211.5 


60.9 


8091.9 




a" 


7835.4 


7793.8 


10081.0 


41.6 


-2287.2 


FbFa 


S 


4698.2 


-552.3 


-595.9 


5250.5 


43.6 




5P 


4759.1 


-491.2 


-535.2 


5250.4 


44.0 




5^ 


-61.0 


-61.1 


-60.7 


0.1 


-0.4 


PbCla 


5 


9624.7 


1356.8 


768.2 


8267.9 


588.6 




(JP 


9732.1 


1464.4 


875.2 


8267.7 


589.2 




5^ 


-107.5 


-107.6 


-107.1 


0.1 


0.5 


PbBra 


5 


12868.7 


2162.5 


1253.4 


10706.2 


909.1 




5P 


13064.2 


2360.7 


1483.5 


10703.5 


877.2 




5^ 


-195.5 


-198.2 


-230.1 


2.7 


31.9 


Pbl2 


5 


18747.1 


3400.5 


2105.7 


15346.6 


1294.8 




<5P 


19022.0 


3677.6 


2439.9 


15344.4 


1237.7 




5^ 


-274.8 


-277.1 


-334.2 


2.3 


57.1 



TABLE VL Isotropic (Kiso) and anisotropic (AK) reduced 
one- and two-bonds spin-spin coupling constants of the lead 
halides and their scalar and spin-dependent contributions cal- 
culated with nonrelativistic Levy-Leblond (NR) and relativis- 
tic spin-free and Dirac-Coulomb (Rel.) Hamiltonians. All 
values are given in lO^^T'^ • J^^. 





Contributior 


Pb-X 


X-X 


Molecule 


i" 'Ki,o 


AK 


^K,so 


AK 


PbFa 


NR 


-652.02 


504.44 


5.21 


13.86 




scalar 


-2397.07 


-1236.65 


1.01 


10.00 




SO 


-661.17 


2262.58 


5.62 


25.72 




Rel. 


-3710.26 


1530.37 


11.84 


49.58 


PbCl2 


NR 


-718.39 


536.92 


-2.09 


47.94 




scalar 


-2427.97 


-1980.39 


7.02 


10.72 




SO 


-863.33 


-841.05 


14.21 


52.57 




Rel. 


-4009.69 


-2284.52 


19.14 


111.23 


PbBra 


NR 


-1261.10 


-991.00 


-22.74 


218.95 




scalar 


-3384.23 


-3011.64 


23.20 


67.45 




SO 


-1573.86 


-1282.44 


58.80 


165.90 




Rel. 


-6219.19 


-5285.08 


59.26 


452.30 


Pbia 


NR 


-1778.76 


-1521.95 


-83.71 


375.45 




scalar 


-4101.85 


-6170.91 


26.72 


326.95 




SO 


-2042.77 


-2348.63 


176.05 


218.92 




Rel. 


-7923.38 


-10041.49 


119.06 


921.32 
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FIG. 2: Relativistic one- and two-bonds reduced isotropic (-R'iso) and anisotropic (Aif) spin-spin coupling constants as a 
function of the halogen nuclear charge Zx ■ 
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